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We show that star formation in galaxy bulges is self-regulating through momen- 
tum feedback, limiting the stellar bulge mass to cx cr"*. Together with a black hole 
mass Mbh « ct^ set by AGN momentum feedback, this produces a linear Mbh — Mb 
relation. At low redshift this gives Mbh/M^ 10"'^, close to the observed ratio. We 
show that AGN feedback can remove any remaining gas from the bulge and terminate 
star formation once the central black hole reaches the AfsH — cr value, contrary to 
earlier claims. We find a mild upward deviation from the cr^ law at higher redshift and 
at higher a. 
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1 INTRODUCTION 

It is now widely accepted that the centre of every medium- 
to-large galaxy contains a supermassive black hole (SMBH). 
The mass M of this SMBH is observed to correlate strongly 
with properties of the host galaxy in two ways. First is 
the correlation with the stellar velocity dispersion a of the 
galaxy bulge. 



Mbh ==; 1.5 x IOVIooMq 



(1) 



(cf. iFerrarese fc Merritt|[2000l : iGebhardt et al.ll2000l '). where 
o"200 = o"/200 kms~^. Second is the correlation with bulge 
stellar mass Mb 



Mb 



1.6 X 10"'' Mi, 



(2) 



(cf. iHaring fc Rb3 120041 '). The first of these relations © is 
considerably tighter than the second ((^J, but both are highly 
significant statistically. 

The existence of these relations strongly suggests a 
connection between the growth of the SMBH and of its 
host galaxy, despite the huge disparity in their masses. 
The underlying reason is evidently that the accretion en- 
ergy Eh = EaccA/c^ released in building the SMBH (where 
face ~ 0.1 is the energy efficiency of accretion) greatly ex- 
ceeds the binding energy of the bulge Ef, ~ Mjcr^ by a factor 



Eh 
Ei 



360 



(3) 



Thus absorption of even a small fraction of the accretion 
energy released b y the SMBH must significantly affect the 
bulge (lKing|l2003l ). This suggests that the process connect- 
ing SMBH growth and galaxy growth is some kind of feed- 
back from the SMBH. Because galaxies are generally opti- 
cally thin this feedback must be mechanical, and so presum- 
ably involves matter expelled during the accretion process. 
This in turn strongly suggests that feedback is a consequence 
of super-Edd ington accretion. This is plausible because the 
ISoltan' (1982) relation linking the growth of SMBHs to lumi- 
nous accretion, and the low observed fraction of ac tive galax- 
ies among all galaxies (e.g. iHeckman et al.ll2004 ). together 
require that most SMBH growth occurs in such phases. 

The remaining obstacle to explaining the Mbh — cr re- 
lation IT} is the very inefficient mechanical coupling of ac- 
cretion energy implied by Q. This explains why the pro- 
posal bv lSilk fc Reed (|l998t ). that a significant fraction of a 
SMBH's Eddington energy rate 



47rGA'/BHC 



(4) 



(where k is the Thomson opacity) is deposited as feedback 
into the surrounding galaxy, implies SMBH masses of order 



only ~ 10" ^ - 10 Ma rather than ~ 10** M© for a typical 
galaxy (cf iKingjIioiol ). Instead, two further features of the 
problem imply a theoretical relation very close to the ob- 
served one IT]), 



Mbh = M^ 



(5) 
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where fg ~ 0.16 is t h e bar yon fraction relative to dark 
matter (cf iKingj |2003| . 120051 ). These new features are (a) 
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the observation that SMBH accretion is never highly super- 
Eddington, so that each emitted photon scatters on average 
only once, implying that the wind from the SMBH region 
carries momentum at the Eddington rate I/Edd/c; and (b) 
in impacting the interstellar medium of the host bulge, the 
wind shocks and loses almost all its energy through Comp - 
ton cooling by the AGN radiation field. As iKind (|2010l ') 
notes, feature (a) implies the we can treat the problem in 
the single-scattering limit, while feature (b) implies that the 
wind impact creates a momentum-dnven outfiow in the host 

ismQ 

It follows that, for SMBH masses A/bh < Mct, an 
outflow from the vicinity of the SMBH vicinity will remain 
trapped very close to it. Once Mbh reaches the critical 
value Ma, the momentum rate L^dd/c of the wind precisely 
balances the total weight -ifga^/G of the overlying gas 
(which, assuming isothermal mass distributions for both the 
dark matter and baryons, is independ ent of radial extent; cf. 
iNavakshin. Wilkinson, fc Kindlioooh . The equality between 
thrust and weight then gives the theoretical relation ((5]), as 
the outflow propagates through the bulge on its dynamical 
time and presumably halts further accretion on to the 
SMBH. Precisely how this last step occurs is still the most 
unclear p art of any theory of the M bh — o" relation (see, for 
example, INavakshin fc Powerllioiol ). 

Self-regulated SMBH growth may explain the Mbh — o" 
relation, but what of the Mbh — Mb relation ([2])? As noted 
m lKindl|2003l ). Compton cooling becomes inefficient beyond 
a certain radius Rc within the bulge. At a SMBH mass 
Mbh = Mcr the shock rapidly reaches this radius and then 
accelerates to a higher speed because of the thermal expan- 
sion of the shocked wind, which now no longer cools. Thus 
the outfiow changes its character from momentum-driven to 
energy-driven at this point. 

This high-speed outfiow removes any ambient gas from 
the bulge, quenching star formation in the process, and so 
it is tempting to identify the original baryonic mass within 
this r adius as th e bulge mass Mi, , as was argued by iKind 
l|2003t ) and lKind (jioOS)). However, such an identification is 
not entirely satisfactory because the relation of the original 
gas mass to the expected stellar mass of the bulge in this 
model remains unclear - it may be all or nothing. The goal 
of this short paper is to address this deficiency of the model. 

We first recall the results of 

iMcLaughlin. King, fc NavakshinI (|2006l ). who pointed 
out that star formation in the nuclear star cluster (NC) 
of a galaxy may provide enough momentum thrust to 
establish its own A/nc — relation, offset from the SMBH 
Mbh — o" relation. They emphasised that the momentum 
feedback produced by young NCs is actually quite similar 
to that produced by a SMBH accreting at its Eddington 



^ Similar arguments are presented in [Murray et al.1 ||2005| ). who 
obtain a result identical to l(5}. However, they assume that it 
is dust rather than electron scattering that drives the outflow, 
which impUes multiple rather than single photon scatterings. The 
resulting momentum flux scales as r L/c, where r is the optical 
depth, and so a SMBH that is formally sub-Eddington in the 
electron scattering sense can radiate effectively at the SMBH's 
Eddington limit. 



rate but smaller by a factor of ~ 20 for a IChabriej (|2003h 
IMF. This similarity arises because massive young stars 
are observed to radiate near their Eddington limit and to 
produce outfiows at rates almost mat ching that of a bl ack 
hole of the same mass in the model of lKinel l|2003l . l2005l ). 

This paper extends the idea of 

IMcLaughlin. King, fc NavakshinI (|2006f ) by investigating 
momentum feedback from a galaxy's stellar bulge, rather 
than its central NC. In particular, we explore the role stellar 
feedback could play in regulating bulge ma ss Mi, and in 
shaping the Mbh — Mb relation. We note that lMurrav et al] 
(2005) followed a similar approac h in th eir explanation of 
the origin of the lFaber fc JacksonI (|l976l ) relation {L oc a*) 
for elliptical galaxies. There is an important distinction, 
however, between f eedback from a galaxy's bulge and feed - 
back from its NC. iMcLaughhn. King, fc NavakshinI (|2006| ) 
considered instantaneous feedback from young massive 
stars in a NC while they are still on their main sequence. 
In contrast, the shortest timescale in a galaxy bulge is 
the dynamical timescale, w hich is much longer than the 
lifet imes of massive stars (cf. INavakshin. Wilkinson, fc Kind 
12009! ). and so we must treat feedback from the bulge in a 
time-integrated or extended form. As for NCs, the feedback 
produced in this case is proportional to the total mass of 
the bulge - but with a smaller efficiency because the stellar 
population in the bulge is relatively mature. 

We explore these ideas in more detail in the following 
sections. We begin by discussing star formation and the role 
it plays in regulating the masses of the galaxy bulges (§[2]); 
then we examine how gas is swept out of the bulge by black 
hole feedback (§0; and finally, we assess the implications 
of these ideas for the Mbh — Mb relation (§|4}. 



2 STAR FORMATION IN GALAXY BULGES 

We begin our analysis by considering how star formation in 
a galaxy bulge impacts on the mass of virialised gas Mg^vir in 
a dark matter halo. A fraction of the virialised gas mass will 
be converted into bulge stars (A/j,), but once the stellar mass 
in the bulge exceeds a critical value (Mmax), momentum- 
driven feedback from the stars (in the form of winds and 
supernovae) will expel this gas and help to regulate the bulge 
mass by suppressing further star formation. Our goal is to 
establish what mass of stars must form to suppress further 
star formation. As we noted in the introduct i on, a similar 
problem has been explored by 'Murr av et al.l ||200E ) as an 
explanation for the origin of the Fa ber fc Jackson (|l976l ) 
relation for elliptical galaxies, L cx ct'*. 

Massive stars dominate the mass and momentum out- 
fiow &wnstaxd_usters^ 

(cf . i Leitherer. R obert fc Dri ssenI 1 19921 : iLamers fc Cassinellil 
I99C 

For a lSalpeter (195 5) initial mass function (IMF) and 
solar metallicity, the contributions peak at stellar masses 
of 12Mq for supernovae and about 50Mq for momentum 
outfiow via winds. Integrated over time, these two feedback 
mechanisms contribute roughly equally to the momentum 
outfiow from massive star clusters. Consequently, the char- 
acteristic main sequence age for stars dominating momen- 
tum feedback is t^a — 10^ years. This is of order the dynam- 
ical time in inner parts of the dark matter halo and shorter 
than the Salpeter timescale (fs — 4 x 10^ years) on which 
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the SMBH grows. This means that the rate of stellar mo- 
mentum feedback is fixed by the rate at which new massive 
stars for m and replace ones that arc dying or have died (see 
Fig. 12 in lLeitherer. Robert fc Driss en 1992,'). In this regime, 
star formation injects momentum to the ambient medium at 
the global rate 

p. ~ e*cM, (6) 

where e» 10""^ (cf. Leitherer. Robert fc Drisserj|l992l : see 
also §2.2 of lMurrav et al.ll2005l l. 

Hence the total momentum produced by a stellar mass 
M* is simply 



e.cM* 



(7) 



Feedback may inhibit star formation when gas that is not 
locked into stars acquires a typical velocity of order the ve- 
locity dispersion cr. Let Mo be the initial bulge gas mass 
prior to the onset of star formation; we assume that Mo is 
of order the virial value Mg^vir- Feedback inhibits further star 
formation in the bulge when the total momentum injected 
into bulge gas reaches 



Pinh ^ 7?Moo-, 



(8) 



where 77 ~ 1. 77 captures the uncertainty in the amount of 
momentum feedback that the system can absorb before star 
formation is suppressed. For example, if the gas needs to be 
completely expelled from the halo by star formation feed- 
back rather than by SMBH feedback, then we might expect 
r] to be greater than unity because the average escape ve- 
locity from the galaxy will be typically just above o. On 
the other hand, if SMBH expels the remaining gas then 77 
could be less than unity: star formation feedback may sim- 
ply need to slow down collapse of gas into s tars until the 
SMBH grows up t o the A/bh — o" value (e.g. iMurrav et al.l 
I2OO5I : iNavakshin.'W ilkinson. & King 200ii). Observational 
evidence can provide important guidance as to what the 

value of T] might be, as we discuss below. 

The derivation of the Mbh — o" relation in iKind (|2005l ) 
shows that mechanical feedback from accretion onto the 
SMBH is confined to a small region near it initially. This 
suggests that stellar feedback is the most plausible mecha- 
nism by which star for mation can be inhib ited up until this 
point (see also § 5 of lMurrav et al.ll2005l ). Therefore, this 
means that the total stellar bulge mass Mi, cannot exceed a 
value Mniax such that p, = pinh; this means that 



e.cAfinax riMoa 
from (O, which gives 

Mb < Mn,ax ^ vMo — , 
e* c 



(9) 



(10) 



or 



Mb < Mmax ~ 0.677Mocr2oo; (11) 

where (7200 — (T/200kms~^. For the moment we assume that 
t; = 1 but we keep in mind our discussion from above. 

3 SWEEPING OUT THE BULGE 

Inspection of relation (|11[) suggests that a significant gas 
mass Mgas might be present in a galaxy bulge as the stellar 
mass asymptotes towards its limiting value Ms; we would 



expect Mgas to be of order the baryon virial mass Mg,vir. 
The observed Mbh — Mb relation then implies that feedback 
from the central SMBH must be able to sweep away this gas 
as A/bh reaches the value A/o- given by the A/bh — o" relation 
((5|. As the escape speed from an isothermal bulge is ~ 2a, 
the SMBH must supply an energy 



E,sc = iA/gas(2cr)^ 



(12) 



to this gas. 

To check this, we assume that the wind driven by super- 
Eddington accretion on to the SMBH provides thrust given 
by the single-scattering limit, i.e. 

^Edd 



MoutV ' 



(13) 



where Mout is the wind mass outflow rate and v is its 
speed. This limit implies that the wind momentum is similar 
to the original photon momentum, i.e. the scattering optical 
depth is ~ 1 and each photon scatters on average about once 
before escaping (cf. Lam ers fc CassincUi 1999). This is likely 
to hold for SMBH accretion because even the dynamical 
accretion rate is not much larger than the Eddington rate 
Afedd for the hole near i ts M — a mass. That this is so is 
shown explicitly in iKind (|201(]| ). which also shows that 



V ~ eaccC 



(14) 



in these circumstances, where face — 0.1 is the accretion 
efficiency. Accordingly the total energy of the outflow is 

Eont = iAA//BH«' = iAA^BH(eaccC)^ (15) 

where AMbh is the increase in SMBH mass during the 
accretion episode considered. The fraction of Eont supplied 
to the bulge gas now depends on the nature of the shock 
as the wind impacts on it. If the shocked gas cools rapidly 
(fast er than the flow time), th e outflow is momentum-driven 
(e.g. IDvsou fc Williani3ll997h and only a small fraction of 
-Bout is used to drive the outflow. However, this turns out 
to be slightly too small to expel it (as we show in the ap- 
pendix). In contrast, if the shocked gas does not readily cool 
(an energy-driven outflow), the thermal expansion of the 
shocked gas means that almost all of Eout is transferred to 
the gas to expel it. 

IKind (|2003l . 120051 ) has shown that the outflow is 
momentum-driven when very close to the SMBH, but 
energy-driven outside a radius ~ 1 kpc for a typical galaxy 
bulgfl This is because the efficiency with which the quasar 
radiation field Compton cools the shock p receeding the out- 
flow declines with increasing radius (cf. ICiotti fc Ostrikerl 
Il997h . If we require that the thrust of the momentum-driven 
outflow must exceed the weight of the overlying gas, as must 
be the case if the outflow is to escape from the immediate 
vicinity of the black hole, then it is straightforward to show 

^ The precise value of this radius scales as (v/c)^. This can be 
seen by comparing the flow time of the outflow to the Compton 
cooling time in the shock precee ding the out flow, given by equa- 
tions (9) and (8) respectively of iKintJ ||2003|) . with equation (14) 
for the flow velocity corrected by a fact o r of 2 to be consistent 
with the value of A^bh derived in iKind l l2005h . In our analysis 
we assume a wind velocity of ~ 0.1c, which is consis tent with 
observational data (cf. Figure 8 of lTombesi et al]|2010l) . 
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that a AfaH — a of the form given by equation ([5]) must 
follow. 

However the total gas mass that the outflow can ulti- 
mately remove is given by considering the total energy iJout, 
as the outflow rapidly becomes energy-driven. By equating 
-Eout and Ecsc we find that 



AMb 



i—y 

VeaccC/ 



= f .8 X fO-Vioo 



(16) 



i.e. only a small increase in SMBH mass is needed to 
sweep the bulge clear of any remaining gas. We quantify 
this further below. 



This implies that A/bh should tend to increase above 
A/o- with redshift and with a. With A/bh = + AA/bh 
the black hole - bulge mass relation finally becomes 



Mb 



7.3 X 10"^ h{z) 



Mb n 

Assuming h = 0.7 a,t z 



1 + 



0.41a-200 
h{z) 

0, this implies that 



Mb 



Mi 



1.15 X W'ri' 



(22) 



(23) 



for a galaxy with (J200 = 1, close to the observed relation 
([2}. This relation should curve slightly upwards for larger a. 



THE BLACK HOLE 
RELATION 



BULGE MASS 



We have argued that momentum-driven outflows from star 
formation tend to produce a bulge stellar mass 



Mb < 0.6?7 Aig,virCr200- 



(17) 



Here we assume that Mo in equation (|ll|l is of order 
the virialised gas mass within the dark matter halo, i.e. 
A^g,vir ~ fg Mv , where Af„ is the total virial mass, and rj cap- 
tures the uncertainty in the amount of momentum feedback 
the system can absorb before star formation is suppressed. 
For the sake of simplicity, we also assume that both the gas 
and dark matter follow isothermal mass distributions, and 
that the velocity dispersion a of the bulge and its underly- 
ing dark matter halo are the same. Neither of these assump- 
tions are likely hold in detail, but their effect is quantitative 
rather than qualitative. Assuming isothermality, how does 
the virial gas mass Afg,vir vary with velocity dispersion a? 

We make the standard assumption that matter is viri- 
alised within a radius such that the mean density is 200 
times the critical value; for the particular case of an isother- 
mal sphere, this gives 



Rv 



5V2H 



(18) 



where H{z) = 100 /!.(z)kms"^Mpc"^ is the Hubble pa- 
rameter at the given redshift z, with h{z) the dimensionless 
Hubble parameter. Because My = 2a^Rv/G we flnd that 



Mb < 5.1 X 10" r? 



1^/9 0^200 



Mr. 



0.16 /i(z) ®' 
when combined with equation ((5)1, we find that 
M^ ^ 7.3 X IQ-* h{z) 



(19) 



(20) 



in reasonable agreement with, for example, 
iHaring fc Ril (|2004l '). 

We expect deviations from this relation and the simple 
AfBH ~ cr* law, as A/bh may grow above M^ by an amount 
< AA/bh in expelling any remaining gas. Combining (|20p 
with (|16|) - where we assume that Afgas is equal to A/g,vir - 
and (|17p gives the predicted maximum deviation from the 
theoretical relation 



AA/bh < r, 4-, '^200 



(21) 



5 DISCUSSION 

We have assumed that star formation in galaxy bulges is 
self-regulating through momentum feedback acting on gas 
in the galaxy's dark matter halo. This puts a limit on the 
total mass of the stars that can form within the halo (equa- 
tion [T71). Crucially, this maximum bulge mass scales as ct"* 
with the galaxy velocity di spersion. C i ven that the SMBH 
mass also scales as (cf. lKinell2003l . l2005l : iMurrav et al.l 
l2005l ). this results in a linear A/bh — Mb relation, close to 
the observed one (|Haring fc Rixl 120041 ) . This relation arises 
because both the SMBH and the bulge are limited by essen- 
tially the same quantity - the maximum momentum thrust 
that the system can take before the gas is bl own away. 

The best match to the observed relation (|Haring fc Rixl 

l2004h is obtained for the dimensionless parameter rj ~ 0.7. 
In the context of our order of magnitude derivation in i|2] a 
value of rj less than one implies that the bulge is not entirely 
self-sufficient in limiting its own growth, and that a SMBH 
or a NC are required to terminate bulge growth completely. 

Finally, we argued in !j4] that if the SMBH growth 
timescale (few x Salpeter) is shorter than the star forma- 
tion (dynamical) timescale, the growth of the central black 
hole is able to shut off star formation in the bulge early. In 
this case the extra SMBH growth introduces a mild upward 
trend in the A/bh oc relation. 
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APPENDIX 

If one assumes that the outflow is momentum-driven at all 
radii (i.e. efflcient shock cooling everywhere), the outflow 
would exert a thrust given precisely by the weight 4fga*/G 
of the overlying gas of mass Mgas. Then the total work done 
by this thrust on the gas, assumed to extend to some radius 
R, is 



because the gas is in equilibrium with the isothermal dark 
matter halo. Therefore a purely momentum-driven outflow 
would fail by a factor 2 to remove the overlying g as, w hich 
is essentially the result claimed bv lSilk fc Nusseil (|2010l ). 

However, cooling is ineffective once the shock radius is 
> 1 kpc for a typical galaxy and so the outflow becomes 
energy-driven. This is because the shock is Compton cooled 
while it is close to the SMBH and the radius at which the 
Compton cooling time begins to exceed the shock flow time 
varies as 0.75(j|oo kpc, where where we've assumed a wind 
velocity of 0.1c. If the shock is to cool radiatively, it must 
do so by means of thermal bremsstrahlung because of the 
magnitude of the shock temperature (T ~ 10^^ {v/c)^K ~ 
10^° K for V — 0.1c). However, the free- free cooling time is 
long in this regime - ~ 5 x 10*n~^ years where n is the 
electron number density - and it exceeds the shock flow 
time at all interesting radii. For this reason, once the shell 
becomes energy-driven, it remains energy driven. 



G 




'CSC 



(24) 



